Disk and elliptical galaxies within renormalization group improved gravity 
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The paper is about possible effects of infrared quantum contributions to General Relativity on disk and ellipti- 
cal galaxies. The Renormalization Group corrected General Relativity (RGGR model) is used to parametrize 
these quantum effects. The new RGGR results presented here concern the elliptical galaxy NGC 4374 and 
the dwarf disk galaxy DDO 47. Using the effective approach to Quantum Field Theory in curved background, 
one can argue that the proper RG energy scale, in the weak field limit, should be related to the Newtonian 
potential. In the context of galaxies, this led to a remarkably small variation of the gravitational coupling G, 
while also capable of generating galaxy rotation and dispersion curves of similar quality to the the best dark 
matter profiles (i.e., the profiles that have a core)F] 



I. INTRODUCTION 

In RefsP^ we presented new results on the applica- 
tion of renormalization group (RG) corrections to Gen- 
eral Relativity in the astrophysical domain, in particu- 
lar on a possible relation between RG large scale effects 
and dark matter-like effects in galaxies. The resulting 
phenomenological model was named RGGR. These de- 
velopments were directly based on the RG application to 
gravity of RefP, and are consistent with the phenomeno- 
logical consequences of diverse approaches to the subject, 
including the related to the asymptotic safety scenario of 
quantum gravity, see in particular RefsP^. 

Currently, in the context of quantum field theory in 
curved space time, it is impossible to construct a formal 
proof that the coupling parameter G is a running param- 
eter in the infrared. However, this possibility can not 
be ruled out. The possibility of General Relativity being 
modified in the far infrared due to the renormalization 
group (RG) has been considered in different contexts, for 
instancep^"^. The previous attempts to apply this pic- 
ture to galax ies have considered for simplicity point-like 
galaxies (e.g.J^J). We extended previous considerations 
by identifying a proper renormalization group energy 
scale /i and by evaluating the consequences considering 
the observational data of dislP' and ellipticaP galaxies. 
We proposed in RefP the existence of a relation between 
fi and the local value of the Newtonian potential (this re- 
lation was reinforced afterwardi^). With this choice, the 
renormalization group-based approach (RGGR) was ca- 
pable to mimic dark matter effects with great precision. 
Also, it is remarkable that this picture induces a very 
small variation on the gravitational coupling parameter 
G, namely a variation of about 10 -7 of its value across a 
galaxy (depending on the matter distribution). We call 
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our model RGGR, in reference to Renormalization Group 
corrected General Relativity. 

Here we present a brief review of the RGGR achieve- 
ments in galaxies and present new results on the galaxies 
NGC 4374 and DDO 47. The first is a giant elliptical 
galaxy that was first analyzed in Ref P in the context of 
RGGR and MOND, the second is a well known dwarf 
disk galaxy whose results within RGGR are for the first 
time here presented, and it is part of a larger work yet 
to be published^. It should be pointed that disk and 
elliptical galaxies behave as stationary systems that are 
stable due to different physical reasons, disk galaxies are 
essentially axisymmetric bodies supported by rotation, 
while elliptical galaxies are about spherically symmet- 
ric and mainly supported by velocity dispersions. Hence 
models for galaxy kinematics may in principle succeed 
in, say, disk galaxies but fail at the ellipticals ones. It is 
remarkable that the RGGR model is doing fine in both 
cases. 



II. A BRIEF REVIEW ON RGGR 

The gravitational coupling parameter G may behave as 
a true constant in the far IR limit, leading to standard 
General Relativity in such limit. Nevertheless, in the 
context of QFT in cuved spa ce time, there is no proof on 
that. According to RefsP^l, a certain logarithmic run- 
ning of G is a direct consequence of covariance and must 
hold in all loop orders. Hence the situation is as follows: 
either there is no new gravitational effect induced by the 
renormalization group in the far infrared, or there are 
such deviations and the gravitational coupling runs as 

fc_ l3 ,_ = 2,-^ = 2*3- (1) 

Equation ([I]) leads to the logarithmically varying G(^) 
function, 
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where Mo is a reference scale introduced such that 
G(no) = Go. The constant Go is the gravitational con- 
stant as measured in the Solar System (actually, there is 
no need to be very precise on where G assumes the value 
of Go, due to the smallness of the variation of G). The di- 
mensionless constant v is a phenomenological parameter 
which depends on the details of the quantum theory lead- 
ing to eq. ([2]). Since we have no means to compute the 
latter from first principles, its value should be fixed from 
observations. Even a small v of about ~ 10~ 7 can lead 
to observational consequences at galactic scales. Note 
that the first possibility, namely of no new gravitational 
effects in the far infrared, corresponds to v = 0. 

The action for this model is simply the Einstein-Hilbert 
one in which G appears inside the integral, namely,^ 



^RGGR [^?] — 
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In the above, G should be understood as an external 
scalar field that satisfies |2]). Since for the problem of 
galaxy rotation curves the cosmological constant effects 
are negligible, we have not written the A term above. 
Nevertheless, for a complete cosmological picture, A is 
nec essa ry and it also runs covariantly with the RG flow 
of « 

There is a simple procedure to map the solutions from 
the Einstein equations with the gravitational constant Go 
into RGGR solutions. In this review, we will proceed to 
find RGGR solutions via a conformal transformation of 
the Einstein-Hilbert action, and to this end first we write 



G = G + SG, 



(4) 



and we assume 5G/Go -C 1, which will be justified latter. 
Introducing the conformally related metric 
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the RGGR action can be written as 

S RGGU [g}=S EH {g}+0(SG 2 



(5) 



(6) 



where S EH is the Einstein-Hilbert action with Go as 
the gravitational constant. The above suggest that the 
RGGR solutions can be generated from the Einstein 
equations solutions via the conformal transformation ([5]) . 
Indeed, within a good approximation, one can check that 
this relation persists when comparing the RGGR equa- 
tions of motion to the Einstein equations even in the 
presence of matted. 

In the context of galaxy kinematics, standard General 
Relativity gives essentially the same predictions of New- 
tonian gravity. In the weak field limit and for velocities 
much lower than that of light, the gravitational dynam- 
ics can be derived from the Newtonian potential, which 
is related to the metric by 



.9oo = - 1 
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Hence, using eq. ([5]), the effective RGGR potential $ is 
given by 



$ = $ 
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An equivalent result can also be found from the geodesies 
of a test particle^. For weak gravitational fields $ n /c 2 <C 
I (with $jv = at spatial infinity), hence even if 
SG/Gq <C 1 eq. ^ can lead to a significant departure 
from Newtonian gravity. 

In order to derive a test particle acceleration, we have 
to specify the proper energy scale fi for the problem set- 
ting in question, which is a time-independent gravita- 
tional phenomena in the weak field limit. This is a recent 
area of exploration of the renormalization group applica- 
tion, where the usual procedures for high energy scat- 
tering of particles cannot be applied straightforwardly. 
Previously tcP the selection of fj, cx 1/r, where r is 
the dista nce from a massive point, was repeatedly used, 
e g) 6 | 8 | io||i2 | This identification adds a constant veloc- 
ity proportional to v to any rotation curve. Although 
it was pointed as an advantage due to the generation 
of "flat rotation curves" for galaxies, it introduced diffi- 
culties with the Tully-Fisher law^, the Newtonian limit, 
and the behavior of the galaxy rotation curve close to the 
galactic center, since there the behavior is closer to the 
expected one without dark matter, we introduced 
a y, identification that seems better justified both from 
the theoretical and observational points of view. The 
characteristic weak-field gravitational energy scale does 
not comes from the geometric scaling 1/r, but should be 
found from the Newtonian potential &n , the latter is the 
field that characterizes gravity in such limit. Therefore, 
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If / would be a complicated function with dependence 
on diverse constants, that would lead to a theory with 
small (or null) prediction power. The simplest assump- 
tion, /i cx <I>jv, leads to /i cx 1/r in the large r limit; which 
is unsatisfactory on observational grounds (bad Newto- 
nian limit and correspondence to the Tully-Fisher law). 
One way to recover the Newtonian limit is to impose a 
suitable cut-off, but this rough procedure does not solves 
the Tully-Fisher issues^. Another one is to use^ 
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(10) 



where $0 and ot are constants. Apart from the condition 
$0 < 0, in order to guarantee 8G/G0 <C I, the precise 
value of $0 is largely irrelevant for the dynamics, since 
$'(r) does not depends on <f> - The relevant parame- 
ter is a, which will be commented below. The above 



energy scale setting ( 10 1 was recently re-obtained from 
a more fundamental perspective^, where a renormaliza- 
tion group scale-setting formalism is employed. 
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The parameter a is a phenomenological parameter that 
needs to depend on the mass of the system, and it must 
go to zero when the mass of the system goes to zero. This 
is necessary to have a good Newtonian limit. From the 
Tully-Fisher law, it is expected to increase monotonically 
with the increase of the mass of disk galaxies. In a recent 
paper, an upper bound on va in the Solar System was 
derivecP. In galaxy systems, ^a| Galaxy ~ 10 -7 , while for 
the Solar System, whose mass is about 10 -10 of that of a 
galaxy, ^a| Solar system 10 -17 . It shows that a linear in- 
crease on a with the mass (ignoring possible dependences 
on the mass distribution) is sufficient to satisfy both the 
current upper bound from the Solar System and the re- 
sults from galaxies. Actually, in RefPit is shown that a 
close-to-linear dependence on the mass can also be found 
for elliptical galaxies by using the fundamental plane. 

Once the /i identification is set, it is straightforward to 
find the rotation velocity for a static gravitational system 
sustained by its centripetal acceleration^, 



V 2 « V 2 

V RGGR ~ V N 



N 



(11) 



Contrary to Newtonian gravity, the value of the Newto- 
nian potential at a given point does play a significant role 
in this approach. This sounds odd from the perspective 
of Newtonian gravity, but this is not so from the Gen- 
eral Relativity viewpoint, since the latter has no free zero 
point of energy. In particular, the Schwarzschild solution 
is not invariant under a constant shift of the potential. 



Equation (11) was essential for the derivation of the 
RGGR galaxy rotation curves. Since elliptical galaxies 
are mainly supported by velocity dispersions (VD), the 
main equation for galaxy kinematics in this case will not 



be eq.( 11 ), but the following expression for the projected 
(line-of-sight) VD (see alscP), 
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B(x,a,b) = Jl t° -1 (l — t b ~ 1 )dt is the incomplete beta 
function, T is the Gamma function, r stands for the de- 
projected (spherical) radius, R for the projected (line of 
sight) radius, I(R) is the luminosity intensity (its inte- 
gral on an infinity surface gives the total luminosity L of 
the galaxy), £(r) is the luminosity density (found from 
the deprojection of I(R)), f3 is the anisotropy parame- 
ter (assumed constant inside each galaxy, it is zero if the 
galaxy has an isotropic VD profile) and M(r) is the total 
(effective) mass of the system at the radius r. See RefP 
and references therein for additional details. In the case 
of Newtonian gravity without dark matter M(r) would 



stand, within a good approximation, as the internal stel- 
lar mass to radius r, i.e., M«(r). 

For RGGR without dark matter, the total mass inside 
the radius r is also the baryonic mass M*(r). Neverthe- 
less, the non-Newtonian gravitational effects of RGGR 
for spherically symmetric systems can be understood 
from the Newtonian perspective as if the total mass was 
given by the following total effective mass^ 



M(r) = M,{r) + M RGGR (r), 



with 
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III. NGC 4374 AND DDO 47 



NGC 4374 



NGC 4374 is a giant elliptical recently analyzed within 
rggrP. Here we present a variation of the analysis pre- 
sented in that reference, namely we here do not directly 
use the photometric data with the Sersic extension to 
model the stellar mass, instead we here only use the Sersic 
model parameters that best fit the surface brightness of 
this galaxy. See Table [T] and Fig. [TJ This is relevant 
to display that our results presented kP are sufficiently 
robust to small changes on the baryonic mass content. 

In RefP it is shown hat MOND fits better the NGC 
4374 observational data than Newtonian gravity without 
dark matter. However, it is still a poor fit, in particular 
since: i) There is a significant tendency towards a lower 
VD curve at large radii, tendency which is strongly en- 
hanced once the fits consider the expected T + ; ii) if the 
expected T* is not used, the best fit is achieved for tan- 
gential anisotropy with j3 < — 1. Other issues of MOND 
with the giant ellipticals can be found for instance in 
RefP^l. It was also shown that RGGR fit to the data is a 
satisfactory one and outperforms MOND in all the points 
above. 



B. DDO 47 

Here we present part of a new result from a working in 
progress^. The dwarf disk galaxy DDO 47 is cited as a 
paradigmatic galaxy for testing dark matter effects^HH^ 
and in particular its fits using the NFW profile vary from 
bad (using the two NFW halo parameters as free param- 
eters) to a disaster (if one of the parameters is fixed in 
accordance with N-body simulations expectations)^. 

Here we use the same data and conventions used in 
Ref)22l for the baryonic matter and observed rotation 
curve, except for the gas surface density. For the lat- 
ter, we directly use the gas surface density provided by 
Ref™, which leads to a gas rotation curve less smooth 
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NGC 4374 

RGGR without dark matter 

Stellar model W avxlO 7 ft T^/(^^) x 2 xL 

Full Sersic+^oj 15.2 ± 2.4 3.94 ±0.56 20.1 0.96 

Full S6rsie+/3 [ _ lil] 19-0^;° 0.57±?;f 2.3±| | 18.5 0.92 

FuU Sersic±K.IMF±^ [0] 13.5 ± 1.4 4.36 ± 0.27 21.9 0.99 

Full S6rsic+K.IMF+ / 3 [ _ 1] i] 13.7 ±1.9 -O^ol^gg 4.44 ± 0.39 21.2 1.0 



TABLE I. NGC 4374 results for RGGR and assuming that the stellar mass profile is only given by a Sersic profile. This table extends 

a table on this galaxy in Ref.E, see this reference for further details. (1) /3[qj indicates isotropic VD, /9r 1,1] indicates constant anisotropy 

with f3 £ [—1, 1], K.IMF is a reference to Kroupa IMF, and it means that the expected value of Y„ was used (i.e., 4.5 ± I.OMq/ Lq y]Pl 
Note that with this full Sersic profile, for this galaxy, the RGGR fit is slightly better. The changes on the best fit parameters above, in 
regard to those found in Ref!2) are well inside the la uncertainties. Hence, small changes in the baryonic matter are indeed generating 
small changes on the model parameters. 




FIG. 1. Two NGC 4374 mass models with RGGR gravity and without dark matter. The curves refer to mass models composed 
by the stellar component (given by a Sersic profile) and RGGR gravity. The black solid lines are the resulting VD for each 
model, the yellow dashed and blue dotted lines are respectively the stellar Newtonian and non-Newtonian contributions to the 
total VD. One of the models assumes isotropy (/3 = 0), the second assumes f3 G [—1,1]. The vertical dashed line signs the 
radius above which the observational data is considered for the fitting procedure (10 arcsec). 



than the one presented in Refill. The results are not 
significantly sensitive to either of the gas profiles and are 
shown in Fig. [2] 



IV. CONCLUSIONS 

In summary, RGGR is a model based on the theoreti- 
cal possibility that the beta function of the gravitational 
coupling parameter G may not be zero in the far infrared. 
Currently, there is no way to directly deduce this behav- 
ior from first principles, nevertheless eq. ([!]) is a natural, 
if not unique, possibility that has appeared many times 
before in the context of Quantum Field Theory on curved 
space-time. This equation depends on a universal free pa- 
rameter v that can be constrained from experiments and 
observations, with v — corresponding to standard Gen- 



eral Relativity. The eq. ([!]) also depends on an energy 
scale [A, which should be related to the symmetries and 
physical interactions that are being evaluated. Consid- 
ering gravitation effects in stationary systems with weak 
Newtonian potential (<£>7v/c 2 <C 1) and slow particle ve- 
locities (v/c <C 1), it is natural to use a relation of the 
type ([9]). In RefPthe relation (10 1 was first proposed. 



RGGR without dark matter is a model with one phe- 
nomenological free parameter (a) which is capable of 
dealing with the kinematics of diverse galaxies. The a 
relation to other physical parameters we aim to under- 
stand soorP^. 
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FIG. 2. DDO 47 rotation curve fits. The red dots and its error bars are the rotation curve observational data, the gray ones 
close to the abscissa are the residues of the fit. The solid black line for each model is its best fit rotation curve, the dashed 
yellow curves are the stellar rotation curves, the dotted purple curve is the gas rotation curve, the dot-dashed green curve 
(shown in MOND and RGGR plots) is the resulting Newtonian, with no dark matter, rotation curve, and the blue dot-dashed 
curve present in the NFW plots is the dark matter halo contribution to the total rotation curve. Note that both of the NFW 
fits favor zero mass-to-light ratios, which is clearly wrong. We stress that the problems with the NFW halo fitting were not 
found here for the first time (we only found them again from slightly different conventions), see Ref P and references therein 
for further details. The result for RGGR above is a new one. For MOND, similar considerations appeared in the ArXiv version 
of the RefP. 
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